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RESUME 


Ce  rapport  a  pour  but  de  dgcrire  METRIC  2,  un  module 
mathfimatique  viser  3  la  determination  des  pieces  de  rechange. 
Apr&s  une  discussion  de  sa  structure  g€n€rale,  ses  moyens 
probabilistes  et  d' optimisation  sont  sondes.  Enfin,  la  sequence 
d'etapes  3  la  solution  d'un  probl&me  donne  est  mise  en  evidence. 
Cette  sequence  en  effet  est  celle  suivie  par  le  programme 
informatique  utilise  dans  la  procedure  d1 approvisionnement 
initiale  d'un  equipement  majeur  (Logistics  Management 
Instruction  1630) . 


ABSTRACT 

"y 

This  report  describes  METRIC  2,  a  mathematical  sparing 
model  for  repairables.  After  discussing  its  general  structure, 
the  probability  and  optimization  techniques  are  explored. 

Finally  the  sequence  of  steps  required  to  solve  a  given  problem 
is  discussed.  These  are  basically  the  steps  used  in  the  FORTRAN 
program  of  the  model  currently  used  in  initial  provisioning  for 
capital  equipment  in  accordance  with  the  procedures  of  CF 
Logistics  Management  Instruction  1630. 
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FOREWORD 


The  objective  of  this  memorandum  is  to  centralize  and 
to  better  explain  all  the  mathematical  theory  behind  METRIC  2, 
a  complex  initial  provisioning  model  designed  to  spare  repairables 
for  weapon  systems  upon  procurement.  It  has  been  found  in  the 
past  that  analysts  taking  over  responsibility  for  the  model  have 
to  locate  and  draw  together  several  unrelated  reports  in  order 
to  understand  the  mechanics  of  the  model  to  a  satisfactory  level. 
Often  these  reports  are  papers  whose  abridged  contents  are  aimed 
at  experts  in  the  field.  Thus  this  work  will  provide  a  funda¬ 
mental  contribution  to  the  continuity  of  the  support  that  D  Log  A 
must  provide  for  METRIC  2. 

As  a  secondary  objective,  the  consolidation  of  this 
material  into  one  document  provides  an  important  first  step  in 
a  project  whose  mandate  is  to  modify  METRIC  2.  This  work  will 
supply  a  solid  foundation  from  which  future  changes  to  the 
model,  consistent  with  its  current  functioning,  can  be  made. 

The  documentation  of  the  mathematical  theory  in  this 
memorandum  is  intended  to  be  self-contained.  The  four  principal 
parts  (the  general  description,  the  probability,  the  optimization 
and  the  computations)  were  designed  to  give  a  good  working 
knowledge  of  the  model's  theory  without  proofs  or  discussions 
of  the  more  involved  aspects.  Annexes  A  through  D  describe 
the  more  complex  mathematical  tools  or  proofs  used  in  the  main 
text.  Those  whose  needs  or  interests  require  more  depth  are 
referred  to  the  respective  section.  A  first  course  in  calculus 
and  probability  and  statistics  is  all  that  is  required  as 
background. 
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METRIC  2:  THE  MATHEMATICAL  THEORY 
OF  A  SPARING  MODEL  FOR  REPAIRABLES 

I  -  GENERAL  DESCRIPTION  OF  THE  MODEL 

BACKGROUND  AND  OBJECTIVES 


1.  Inventory  control  generally  deals  with  two  questions: 
when  to  reorder  and  how  much  to  reorder.  However,  there  is  a 
significant  class  of  items  where  demand  is  so  low  or  cost  so 
high  that  reordering  immediately  upon  demand  is  the  best 
possible  policy.  The  Canadian  Forces  stock  contains  such  a 
class  -  the  repairable  or  recoverable  items.  Although  they 
make  up  six  percent  of  the  total  stock,  their  dollar  value  is 
about  45  percent  of  the  total  investment  at  any  given  time. 

In  order  to  satisfy  demands ,  a  spare  stock  of  these  items  should 
be  kept  to  guard  against  stock-outs,  for  resupply  times  are 
non-zero.  Thus  the  optimal  inventory  policy  (spare  stock)  must 
consider  the  tension  between  the  high  cost  of  this  spare  stock 
and  the  risk  of  stock  out.  This  is,  in  broad  terms,  the  nature 
of  the  problem  treated  by  METRIC. 

2.  Before  describing  this  situation  any  further  it  is 
imperative  that  clear  definitions  of  terms  be  made  here.  The 
term  reorder  will  be  taken  to  mean  a  request  for  repair.  An 
item  is  in  resupply  if  it  is  in  the  process  of  being  repaired. 
Let  x  be  the  number  of  units  on  order  at  a  random  point  in  time 
and  let  8  be  the  total  spare  stock  level.  Then  if  x  <  s, 
these  x  units  are  said  to  be  in  routine  resupply.  If  x  >  s, 
then  x  -  s  units  are  said  to  be  on  backorder.  Thus  the  total 
spare  stock  equals  the  stock  on  hand  plus  on  order  minus  back¬ 
orders.  As  we  shall  see,  a  driving  factor  in  establishing  this 
spare  stock  level  s  is  the  average  time  required  in  resupply  or 
mean  time  to  repair  (MTTR) .  In  turn  this  resupply  time  is 


directly  affected  by  the  geographical  hierarchy  of  the  repair 
facilities.  This  partitions  the  repair  process  into  echelons . 

For  example  a  unit  removed  fran  an  equipment  (aircraft  for 
example)  will  be  brought  into  the  repair  shop  at  the  base  out 
of  which  this  equipment  operates.  If  the  unit  cannot  be 
repaired  there,  it  will  be  sent  to  a  DND  depot.  Again  if 
repairs  cannot  be  done  there  it  will  be  sent  to  an  outside 
contractor.  Thus  first  echelon  will  be  the  base,  second 
echelon  will  be  an  in-house  depot  and  third  echelon  a  contractor. 
Another  driving  factor  defining  the  stock  level  s  is  the  mean 
time  between  failures  (MTBF) ,  as  this  establishes  the  demand 
rates  and  thus  the  pressure  on  the  spare  stock. 

3.  METRIC,  which  is  short  for  "Multi-Echelon  Technique 

for  Recoverable  Item  Control",  was  developed  as  a  systems 
approach  to  initial  provisioning  of  repairable  items  for  a 
weapon  system.  In  particular,  the  total  investment  in  repairable 
spares  and  their  distribution  is  related  to  the  efficiency  with 
which  they  support  the  major  equipment  of  the  weapon  system. 

This  permits  life  cycle  costing  and  comparisons  of  different 
alternatives  when  purchasing  a  weapon  system.  In  particular 
with  estimates  of  failure  rates  and  servicing  rates  under  given 
economic  or  operational  constraints,  it  will  optimize  a  given 
measure  of  support  function.  Also,  given  existing  stock  levels 
at  various  bases  and  depots,  it  can  evaluate  the  level  of 
support  provided  by  this  allocation  and  can  redistribute  the 
stock  to  optimize  the  measure  of  support.  METRIC  optimizes 
only  in  terms  of  what  stock  is  bought  and/or  how  it  is 
distributed.  It  does  not  optimize  for  instance  against  repair 
policy,  scheduling  or  transportation  systems  between  bases 
and  depots.  In  what  follows,  the  word  equipment  will  refer  to 
the  major  equipment ,  a  certain  number  of  which  form  the  weapon 
system.  An  equipment  is  assumed  to  be  made  up  of  a  certain 
number  of  items  (type) .  A  particular  repairable  piece  of  a 
certain  item  type  will  be  referred  to  as  a  unit  of  that  item. 
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Also  in  what  follows  bracketed  numbers  (e.g.  (3))  in  the  text 
give  the  pertinent  references.  The  symbol  □  will  refer  to  end 
of  proofs. 

THE  METRIC  MODEL 

4.  The  METRIC  model  has  a  two  echelon  structure,  that 

is  there  are  a  certain  number  of  bases  being  supported  by  one 
central  depot  or  by  a  series  of  depots  each  responsible  for 
servicing  certain  bases  or  for  repairing  certain  items  only. 

When  a  unit  comes  into  a  base  for  repair  it  has  a  certain 
probability  r  of  being  repaired  at  the  base  and  a  probability 
1-r  of  being  repaired  at  a  depot  (see  figure  1) .  The  two 
basic  parameters  of  the  model  determining  the  flow  of  repair- 
ables  through  the  system  are  the  MTBP  and  the  MTTR.  The 
model  assumes  that  demands  are  governed  by  a  compound  Poisson 
distribution  with  a  certain  variance  to  mean  ratio  which  is 
identical  for  all  items  (see  II  -  Demand  Process).  The  MTBF's 
provide  the  customer  arrival  rates  for  the  compound  Poisson 
distribution  and  a  theorem  of  Palm  tells  us  we  can  use  the  MTTR 
as  the  time  period  in  the  compound  Poisson  distribution  in 
order  to  calculate  the  state  probabilities  (that  x  units  are 
in  resupply  at  a  random  point  in  time) .  The  MTTR  is  a  composite 
of  the  base  repair  time  (BRT) ,  the  depot  repair  time  (DRT)  and 
the  order  and  shipping  time  (OST) .  Prom  these  probabilities 
the  measure  of  performance,  availability,  can  be  computed  for 
each  item.  The  availability  of  an  item  is  simply  the  probability 
that  a  demand  can  be  satisfied  by  the  stock  on  hand,  i.e.  the 
probability  that  no  backorders  exist  at  a  given  point  in  time 
(see  II  -  Measure  of  Supply  Performance) .  As  the  failures  of 
the  items  comprising  an  equipment  are  assumed  to  be  independent, 
the  product  of  the  items'  avilabilities  gives  the  probability 
that  no  backorders  exist  on  any  item  at  a  random  point  in  time, 
i.e.  the  fraction  of  time  all  equipments  are  capable  of  operating. 
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This  is  called  the  system  availability.  If  there  are  n  equip¬ 
ments  operating,  the  availability  of  one  equipment  chosen  at 
random  can  be  obtained  by  taking  the  n ^  root  of  the  system 
availability. 


5.  At  initial  provisioning  estimates  of  the  MTBF's  are 
at  best  uncertain.  The  METRIC  model  compensates  for  this  by 
using  a  Gamma  distribution  for  the  true  demand  with  the 
estimate  of  MTBF  as  its  mean.  This  distribution  is  integrated 
with  the  availability  function  to  obtain  a  mathematical  expect¬ 
ation  of  the  latter.  If  demand  data  are  available,  observations 
will  be  used  to  modify  this  gamma  distribution  using  a  Bayesian 
technique  (see  II  -  Demand  Prediction) . 

6.  The  METRIC  model  uses  marginal  analysis  to  obtain  the 
most  efficient  allocation  of  stock  or  investment.  It  can 
operate  in  three  different  modes.  First,  given  an  availability 
target  for  each  equipment  it  will  minimize  the  investment,  or 
given  an  investment  target  it  will  maximize  the  availability 

for  each  equipment.  Table  3  of  (3)  is  reproduced  here  in  Table  I 
and  gives  two  examples  of  METRIC  calculations,  one  with  an 
investment  constraint  and  one  with  an  availability  constraint. 

In  each  it  is  assumed  that  two  equipments  are  located  at  each 
of  the  four  bases  all  with  the  same  operational  activity  per 
year  (see  III).  Second,  it  can  redistribute  an  existing  stock 
distribution  so  as  to  optimize  the  availability.  Third,  it  can 
evaluate  an  existing  stock,  that  is,  it  can  determine  the  level 
of  availability  for  a  given  stock  distribution  (see  IV,  para¬ 
graph  2)  • 
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T/.I'LE  I 

EXAMPLES  OF  METRIC  CALCULATIONS 


MONEY 

CONSTRAINT  OF  $25,000 

BASE* 
REPAIR 
ITEM  FRACTION 

BASE 

REPAIR 

TIME 

(MON) 

DEPOT 

REPAIR 

TIME 

(MON) 

UNIT 

COST 

($) 

ESTIM. 

FAILURES 

PER  YEAR 
(SYS) 

METRIC’S 

STOCKAGE 

POLICY 

ITEM 
.  AVAIL. 

% 

TOTAL  PIPELINE  BASES 
12  3  4 


1 

100% 

.50 

0 

100 

9.63 

8 

0 

2 

2 

2 

2 

.  9964 

2 

50% 

.50 

2.0 

400 

4.89 

6 

2 

1 

1 

1 

1 

.9877 

3 

20% 

1.0 

2.0 

1000 

7.74 

7 

3 

1 

1 

1 

.9728 

4 

0% 

0 

2.0 

3000 

3.00 

S 

1 

1 

1 

1 

1 

.9805 

Total  Investment  -  $25 <200 

System  Availability  -  93.86 

Equipment  Availability  -  (93.87)1/,S  »  99.21% 


EQUIPMENT  AVAILABILITY  CONSTRAINT  OF  95% 


ITEM 

BASE* 

REPAIR 

FRACTION 

BASE 

REPAIR 

TIME 

(MON) 

DEPOT 

REPAIR 

TIME 

(MON) 

UNIT 

COST 

($) 

ESTIM. 
FAILURES 
PER  YEAR 
(SYS) 

t 

METRIC'S 

STOCKAGE 

POLICY 

ITEM 

AVAIL. 

% 

TOTAL 

PIPELINE 

BASES 

1 

2 

3 

4 

1 

100% 

.50 

0 

100 

9.63 

4 

0 

1 

1 

1 

1 

.9647 

2 

50% 

.50 

2.0 

400 

4.89 

5 

1 

1 

1 

1 

1 

.9792 

3 

20% 

1.0 

2.0 

1000 

7.74 

4 

0 

1 

1 

1 

1 

.8244 

4 

0% 

0 

2.0 

3000 

3.00 

1 

1 

0 

0 

0 

0 

.8170 

Total  Investment  -  $9,400 

System  Availability  -  63.62% 

Equipment  Availability  -  <63.62)1/,<I  -  95% 

•(%  of  failed  items  which  can  be  repaired  at  base  level) 
t  The  number  of  applications  of  each  item  in  the  equipment 
and  the  total  annual  operation  hours  are  accounted  for  here. 

H.B.  It  is  assumed  that  the  equipment  and  activity  ere  equally  distributed 
across  the  bases. 
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THE  MODEL'S  PARAMETERS 

7.  Some  of  the  parameters  have  been  mentioned  already. 

We  list  here  the  primary  parameters  required  for  the  model. 


a.  By  Systems 


-  number  of  bases 

-  number  of  items 
number  of  equipments 

change  in  level  of  operation  from  data  period 
to  prediction  period 


b.  By  Item 


-  unit  cost 

-  MTBF 

-  redundancy  within  an  equipment 

-  observed  demand 


number  of  months  of  observed  demand 


c.  By  Item  and  Base  or  Depot 

-  average  base  repair  time  (BRT) 
average  depot  repair  time  (DRT) 

fraction  of  failures  that  are  base  repairable 
average  ordering  and  shipping  time  between 
a  base  and  a  depot  (OST) 

-  number  of  equipments  at  each  base 

-  minimum  and/or  maximum  stock  levels. 


•nrm' 
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THE  MODEL'S  ASSUMPTIONS 


8.  Next  we  list  the  major  structural  assumptions  of 

the  METRIC  model. 

a.  Demand  Distribution:  demand  for  each  item  is 
assumed  to  be  a  logarithmic  Poisson  process 
with  a  variance  to  mean  ratio  assumed  constant 
across  all  items.  Demand  is  also  assumed 
stationary  over  the  prediction  period  though  it 
is  possible  to  vary  the  level  of  operations 
from  the  data  period. 

b.  Repair  Decision;  the  decision  to  repair  a  unit 
at  base  level  or  depot  is  a  function  only  of 
the  type  of  malfunction  and  base  maintenance 
capability.  However,  because  of  theoretical 
constraints  it  is  assumed  that  a  customer  will 
have  all  his  demands  repaired  at  base  level  or 
all  at  the  depot.  Finally,  at  the  depot  it  is 
assumed  that  service  is  done  on  a  first-come 
first-served  basis  and  that  no  batching  of  items 
for  repair  is  permitted. 

c.  Lateral  Resupply:  lateral  resupply  between 
bases  is  not  taken  into  account  in  the  model. 
Resupply  at  a  given  base  in  the  model  comes  from 
the  base  maintenance  shop  if  repair  is  done  at 
that  base;  otherwise  it  comes  from  the  depot. 

d.  No  Condemnations;  the  condemnation  rate  is 
assumed  to  be  zero;  that  is,  an  item  is 
repairable  at  base  or  depot  regardless  of  its 
age  or  previous  number  of  repairs. 
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e.  All  items  comprising  an  equipment  are 

considered  equally  important  to  the  operation 
of  the  equipment  and  are  independent  from  one 
another  (probabilistically  in  their  failures). 

REFERENCES 


9.  The  original  METRIC  model  was  developed  by 

C.  Sherbrooke  (28).  Shortly  afterwards,  a  computer  program 
was  put  together  (8) .  A  non-technical  description  can  be 
found  in  (29) .  This  model  as  developed  by  Sherbrooke  used 
expected  number  of  backorders  as  its  performance  measure. 

When  the  CF  inherited  the  model,  preference  was  indicated  for 
availability  as  the  measure  of  performance.  Studies  have  been 
done  (31)  ,  (20)  indicating  that  the  difference  in  stock 
allocations  incurred  by  switching  from  the  "backorder"  to  the 
"availability”  measure  is  marginal.  The  changes  were  made 
(13)  and  the  version  resulting  was  called  METRIC  2.  During 
the  period  of  implementation  several  reports  were  produced 
on  applications  to  actual  purchases  (  (6) ,  (15)  ,  (16)  ,  (19)  ) 
and  various  studies  were  performed  on  the  validity  of  the  model 
(  (3),  (4),  (5),  (9)).  Though  it  has  been  difficult  to  validate 
the  model  and  though  it  is  clear  that  there  are  shortcomings, 
there  are  at  present  no  alternative  means  that  take  a  systems 
approach  to  the  problem  of  sparing  in  the  acquisition  of 
capital  equipment.  In  response  to  these  shortcomings,  work 
has  been  done  to  get  METRIC  to  behave  more  closely  to  the 
environment  in  the  Canadian  Forces  Supply  System  (CF SS) . 
Adaptation  for  more  bases,  more  spares  and  more  echelons  were 
considered  (  (18)  ,  (22)  ) .  Modifications  to  include  Indent 
levels  (or  part  hierarchies)  have  been  studied  as  well  (  (21) , 
(1)  ) .  At  the  moment  METRIC  is  used  exclusively  for  initial 
provisioning  of  spares.  There  is  a  need  for  METRIC'S 
capability  to  follow  sparing  into  the  reprovisioning  phase 


10 


of  the  CFSS,  as  spares  bought  initially  on  METRIC's 
allocation  should  be  re-evaluated  on  the  same  basis.  Some 
initial  consideration  has  been  given  to  this  problem  (2) . 
Recently  work  outside  DND  has  been  done  on  problems  of  indent 
level,  three  echelons,  finite  repair  capacity  or  repair 
facility  breakdown  (  (24) ,  (14) ,  (7) ,  (17) ,  (23) ,  (25)  ) . 

II  -  METRIC'S  PROBABILITY  MODEL 


DEMAND  PROCESS 


10.  Customer  demand  for  each  item  is  assumed  to  be  a 

logarithmic  Poisson  process.  This  is  a  member  of  the  compound 
Poisson  family.  These  terms  will  now  be  defined.  The  follow¬ 
ing  discussion  refers  to  the  demand  process  of  one  item  at 
base  level.  Capital  letters  will  be  used  to  denote  random 
variables  and  the  corresponding  lower  case  letters  will  re¬ 
present  values  taken  on  by  the  random  variable. 


11.  In  a  Poisson  process  the  distribution  of  customer 

arrivals  is  given  by  the  following  probability  function: 


P(xixt)  =  »t)X  «'*! 

x! 


This  is  the  probability  that  x  customers  arrive  in  an  arbitrary 
time  interval  of  length  t.  This  formula  can  be  derived  from 
the  following  two  assumptions: 


a.  the  probability  of  a  single  arrival  in  a 
small  time  At  is  XAt  (X  =  "arrival  rate") j 

b.  the  probability  of  more  than  one  arrival  in 
At  is  negligible. 
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Note  that  t  is  not  absolute  time  but  rather  the  length  of  an 
arbitrary  interval  of  time.  Here  each  customer  arrives  with 
one  demand  and  hence  the  demand  process  is  said  to  be  a  Poisson 
demand  process.  The  mean  of  this  distribution  is  Xt  and  this 
value  is  also  its  variance.  Finally  the  mean  time  between 
arrivals  (MTBA)  (or  mean  time  between  failures  (MTBF)  as 
failures  generate  arrivals)  is  i  .  These  are  standard  results 
which  can  be  found  in  text  books  dealing  with  this  subject 
(see  (27;  sec.  2-4)  for  instance). 

12.  In  a  compound  Poisson  process  customers  arrive  as 

in  a  Poisson  process  at  a  rate  X  with  one  or  more  demands 
given  by  a  discrete  positive  random  variable  W  with  probability 
distribution  function  f (w)  called  the  compounding  distribution 
function.  Thus  in  an  interval  of  length  t,  the  total  number 
of  demands  is  given  by; 

N  ( t) 

x(t)  =  E  w  i 

i*=l  1 

where  N(t)  is  the  number  of  customers  arriving  in  that  interval, 

as  given  by  the  Poisson  process  and  {wj.  .  are 

1  f  •  ;N(t) 

independent  identically  distributed  random  variables  of  the 
demands  of  customers  arriving  in  the  interval.  We  would  again 
like  to  describe  the  probability  P(x|Xt)  =  Pr(X(t)=x)  of  x 
demands  occurring  in  an  interval  of  length  t  when  customers 
arrive  at  a  rate  X.  This  time,  since  each  customer  must  have 
at  least  one  demand  the  probability  of  y  customers  arriving 
with  x  demands  (y  4  x)  is 

Pr (N (t) *=y  and  X(t)=x)  =  Pr(N(t)*=y)  •  Pr(X(t)*x|N(t)«=y) 

=  (Xt)y  e~Xt  .  Pr(  Z  W,  = 
y!  i=l  1 


x) 


y 

In  the  last  term  we  note  that  for  a  fixed  y,  the  variable  I  V 
has  a  distribution  called  the  y-fold  convolution  of  f  ^  1 
(the  distribution  of  W)  with  itself.  The  probability  that  the 
variable  takes  the  value  x  is  denoted  fy  (x) .  Thus  to  form 
P(x|Xt)  we  must  add  these  probabilities  over  all  y's 
producing  x  demands,  i.e. 


P (x I Xt) 


y* 

f  (x) 


To  see  how  this  distribution  varies  from  the  Poisson  we  can 

2 

calculate  the  mean  m  and  variance  a  (for  details  see  Annex  A) 
They  are: 

m  =  Xt  E (W) 

a2  =  Xt  E(W2) 
o 

where  E(W)  and  E(W  )  are  the  first  two  moments  of  f  given  by 


E(W)  =  Z  w  f(w) 
w=l 


E(W2)  ■  I  w2  f(w) 


As  the  difference  between  the  Poisson  and  the  compound  Poisson 
is  caused  by  f  a  better  measure  of  the  difference  would  be  to 
take  the  variance  to  mean  ratio  q  as  it  eliminates  the  parameter 
Xt  of  the  Poisson  leaving  us  with  the  following: 


q  “  5“ 


-  w-1  v 

to 

I  w  f  (w) 

W"1 


Note  that  we  have  q  *  1  if  f(w) 
the  compound  Poisson  is  Poisson. 


«=  0  for  all  w  >  2,  that  is  if 


-  13  - 


13.  We  must  now  specify  a  distribution  f (w) .  The 

following  function  is  used  as  it  gives  rise  to  a  convenient 
computational  form  for  P(x|Xt)  and  when  q<  3  remains  faithful 
to  previously  accepted  distributions  (q  will  be  shown  to  be 
the  same  as  the  q  in  paragraph  12  above) : 

10  w  =  0 

-4—  (A  w 

w In  q  \q ) 

We  specify  also  that 

U  =  k  In  q  k  >  0 


W  a  1,2,3/  ... 
q  *  p  +  1  >  1 


With  f(w)  defined  as  such,  the  compound  Poisson  process  is 
called  a  Logarithmic  Poisson  Process.  From  the  Taylor  series 
expansion  for  In  we  have: 

*  -fcn  (1  -  -)  =  £,n  q 

w=l  "  q 


so  that  E  f(w)  =  1  as  required  by  a  p.d.f.  Also  the  second 
w=0 

defining  equation  above  gives  rather  simple  expressions  for 
the  mean  and  the  variance  to  mean  ratio  of  the  compound  Poisson 


"  "  Xt  E(M)  *  k  ln  13  ^  wHEihf 


1 


\ 


_  E 
q 


q 


i) 


*» 


-  14  - 


2 

o 

m 


E(W2) 
E  (W) 


E  w-2..-  (e\w 

w=l  w  in  q  \qj 


^3zL 

Inq 


=  Z 
w=l 


Z 


w 


W=1 


w 


P 


The  second  last  equality  results  from  the  fact  that  the 
infinite  sum  is  the  expected  value  of  a  geometric  distribution 
with  probability  of  success  Thus  the  parameter  q  defining 
f (w)  is  the  variance  to  mean  ratio  of  the  final  compound 
Poisson  process.  Through  the  use  of  characteristic  functions 
it  is  shown  in  Annex  A  that  P(x|Xt)  is  a  negative  binomial 
distribution/  namely  if  q  and  k  are  the  parameters  of  the 
logarithmic  Poisson 


where  by  convention 


jr.tt.-i> 


to  permit  y  to  be  any  real  number. 


...  (y-x+1)  is  extended 
x! 


The  calculations  of  P(x)Xt)  can  be  done  recursively 

(x+k)  _  , 

P(x+ljxt)  =  -  £  P(x|Xt), 

(x+1)  * 


k 


where  P(0|Xt) 


MEASURE  OF  SUPPLY  PERFORMANCE 


14.  There  are  several  different  measures  of  supply 

performance  that  can  be  used.  The  original  Rand  version  of 
METRIC  uses  expected  number  of  backorders  at  a  random  point 
in  time.  METRIC  2  uses  the  concept  of  availability  (see  (31)) 
the  probability  of  no  backorders  at  a  random  point  in  time. 
These  values  are  determined  by  the  probability  distribution 
of  the  random  variable  X ,  the  number  of  units  in  resupply. 

The  steady  state  probability  distribution  h(x)  for  X  is  given 
by  a  theorem  due  to  Palm  (see  Annex  B)  which  says  that  if 
demand  is  Poisson  then  X  is  also  Poisson  which  depends  on  the 
mean  of  the  resupply  time  distributions  (mean  time  to  repair  - 
MTTR)  and  not  on  the  resupply  time  distribution  itself.  More 
precise ly : 

Theorem  (Palm)  -  Let  the  demand  for  an  item  be  Poisson  with 
rate  X  and  the  resupply  time  be  an  arbitrary  probability 
distribution^  (t)  with  mean  T.  Then  the  steady  state 
probabilities  of  x  units  in  resupply  are  given  by  the  Poisson 
with  rate  XT,  i.e. 

h(x)  =  e-  "■ T- 

x: 


Thus  the  distribution  h(x)  of  the  number  of  units  in  resupply 
depends  only  on  the  mean  T  of  the  distribution  f(t)  of  the 
resupply  times  and  not,  for  instance,  on  the  variance  of  ^(t). 

15.  To  adapt  this  to  our  situation  of  compound  Poisson 

demands,  we  must  assume  that  a  resupply  time  drawn  from  ^ ( t) 
is  applicable  to  all  demands  placed  by  that  customer.  In  this 
case  then  the  steady  state  probability  that  y  customers  are 
in  the  system  is  given  by 

.  UT)*  «~XT 


h( 
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Now  to  gear  this  to  demands,  since  the  resupply  process  is 
independent  of  the  customer  order  size,  this  must  be 
multiplied  by  the  probability  that  y  customers  place  x  demands, 
i.e.  by  f 1  (x) ,  and  summed  over  all  values  of  y  <  x.  Thus  the 
probability  that  x  units  are  in  resupply  is 

P (x |  XT)  *  Z  (Wy  e  fy*(x) 

y=o  y' 

and  we  get  that  the  number  of  units  in  resupply  is  also 
compound  Poisson  depending  on  T  and  not  ^  (t)  . 

16.  The  availability  is  the  probability  that  no  back¬ 
orders  occur,  i.e.  Pr(X  -4  s)  where  s  is  the  total  spare  stock. 
Denoting  this  by  A(s)  we  have 

s 

A (s)  *  Z  P(xjAT) 
x=0 

If  as  in  the  Rand  version  we  were  concerned  about  the  expected 
number  of  backorders  as  a  measure  the  formula  would  be 

09 

B  ( s)  =  Z  (x-s)  P  (x  |  AT) 

x=s+l 

DEMAND  PREDICTION 

17.  In  order  to  use  the  compound  Poisson  we  must 
determine  the  true  mean  demand  0*AL?fora  period  L  where 
f  is  the  expected  value  of  the  compounding  distribution.  It 
is  assumed  to  be  stationary  but  otherwise  is  unknown.  Because 
inaccuracies  of  point  estimates  for  6  are  magnified  in 
computations  this  method  is  not  acceptable.  Furthermore, 
initial  estimates  of  demand  have  a  tendency  to  be  larger  than 
the  true  mean,  usually  by  a  factor  of  two  (4).  Because  of 
this  inability  to  establish  a  true  mean  demand,  a  Bayesian 


I 
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Procedure  is  used  and  will  now  be  described.  This  procedure 
requires  a  prior  probability  distribution  for  the  values  of 
the  true  mean  demand  which  is  assumed  to  be  a  gamma  distribution 
over  some  unit  period  of  time 

g ( 6)  =  . .  1  ev_1  e  “e/,w  o  <  e  <  » 

WV  r<v)  v'w  >  0 

2  2 

with  mean  m  =  vw  and  variance  a  =  vw  .  The  initial  estimate 
of  0  is  taken  as  this  mean.  The  prior  probability  distribution 
will  eventually  be  meshed  together  with  observed  data  for 
each  item.  However,  data  for  different  items  will  be  observed 
over  different  time  periods.  If  we  are  now  interested  in  the 
true  mean  distribution  over  a  period  r  times  as  long,  then 
we  require  the  distribution  of  X  =  r8.  The  Jacobian  of  this 
transformation  is  ^  so  that  the  distribution  of  X  is 

(|)V"  e-V-l 

_  _ 1 _ _  xv~l  e-X/,rw 

v 

(rw)  T(v) 

This  is  again  conveniently  a  gamma  distribution  with  a  mean  r 
times  the  mean  for  the  unit  period  gamma  distribution  and 
variance  r2  times  the  variance  for  the  unit  period.  Thus 
it  suffices  to  establish  the  true  mean  demand  distribution 
for  some  convenient  time  period. 

18.  Since  there  are  two  parameters  to  be  estimated, 

two  estimation  equations  will  be  needed.  As  mentioned 
earlier  one  is  the  mean 

=  vw 


■  g 


V 

w  r  (v) 
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for  each  item  i.  In  the  absence  of  further  data  we  can  turn 
once  again  to  the  variance  to  mean  ratio.  From  experience 
an  estimate  of  this  value  across  all  items  seems  appropriate 
though  the  above  mentioned  adjustments  for  different  time 
periods  must  be  made.  This  ratio  will  be  called  the  prior 
distribution  ratio  and  the  computer  program  will  be  run  for 
several  different  values  of  the  ratio  for  comparison  of  the 
different  resulting  stockage  policies. 

19.  Once  the  prior  distribution  g(0)  for  an  item  is 

identified  we  want  to  modify  the  distribution  to  take  into 
account  the  fact  that  u  demands  for  the  item  have  been 
observed  in  time  L.  This  is  done  using  Bayes  theorem  to 
obtain  the  following  posterior  distribution  <J>(8|u)  for  0 

.  P(u| 6/f)  g(0) 

<f>(e|u)  *  - 

/p(uU/T)  g(5)  d£ 

C 

Division  by  the  mean  number  of  demands  per  customer  f 

changes  the  mean  demand  rate  to  the  mean  customer  arrival 

rate  as  is  required  by  the  compound  Poisson  probability 

function  P(x|Xt).  Computations  of  availability  A(s)  (see 

Section  IV)  are  done  with  a  time  period  T,  the  average 

response  time.  Thus  the  mean  arrival  rate  0/F  calculated 

T 

on  a  time  period  L  must  be  re-scaled  by  j-  for  the  formula  of 
A (s) .  Furthermore,  if  there  is  a  change  in  the  level  of 
operations  by  a  factor  of  a  from  the  data  period  to  the 
prediction  period,  this  factor  should  also  multiply  the 
customer  arrival  rate.  Thus 

A(s|  8)  =  I  P(x|  a.  1  .  |  )  . 

x**0  L  * 
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Finally  the  values  of  A  ( s  |  0 )  are  weighted  usings  (6|u)  to 
obtain 


A* ( s | u) 


4>(0  |u)  A ( s  !  6) 


de  . 


20.  The  above  integrals  are  approximated  using  sums 

on  five  to  ten  points.  Each  subinterval  is  referred  to  as  a 
Bayes  state  and  in  the  nfc^  Bayes  state  an  appropriate  value 
0n  is  chosen.  This  leads  to  the  following  approximations 

P(u| 0  /f  )  g(0  ) 
u)  —  rc  n 

Z  P(u|0k/f)  g(0k) 

1c 


♦  (enl 


and 


s 

A* (s|u)  «  Z  <j>(elu)  Z 

n  n  x=o 


P(xja 


T  6n 

'n'f 


) 


Note  that  <p  and  g  are  no  longer  probability  density  functions 
but  the  derived  state  probability  obtained  from  suitable 
partitioning  of  the  0  axis.  In  the  absence  of  demand  data 
<H0n|u)  reduces  to  g ( ©n)  .  The  evaluation  of  P(x|X)  requires 
the  two  parameters  q  and  k.  The  value  of  q  is  chosen 
usually  between  1  and  2.  Experience  indicates  (4)  that  a 
value  of  1.1  or  1.2  is  appropriate  in  the  CF.  As  the  mean 
of  P ( x | A )  is  k(q-l)  then  in  each  Bayes  state  division  of  6n 

by  q-1  will  yield  the  appropriate  value  of  k  for  that  Bayes 
state. 


20 


III  -  THE  OPTIMIZATION  TECHNIQUE 

THE  OBJECTIVE  FUNCTION 

21.  Let  =  unit  cost  of  item  i,  i=l,2, - ,  I  where 

I  is  the  number  of  different  items 

a  =  minimum  acceptable  system  availability 
target 

siQ  =  depot  stock  for  item  i 

s^j  =  base  j  stock  for  item  i,  j=l,2, - ,J 

where  J  is  the  number  of  bases 

A(siQ  ,  s^)  =  availability  of  item  i  at  base  j  when  base 
stock  is  sij  and  depot  stock  is  s^Q. 

The  objective  is  to  solve  the  following 
I  J 

min  E  c  Z  s.. 

{si;.}  i=l  j=0  3 

subject  to 

J  I 

n  n  A(sio  '  sijJ  *  a 

j=l  i=l 

or  more  conveniently 
J  I 

E  E  log  A ( 8 .  n  ,  8.  .)  £  log  a 
j=l  i=l  10  13 


If  budgetary  constraints  are  of  primary  concern  the  problem 
will  be  formalized  as 


J  I 

max  I  Z  log  A ( s . -  ,  s . . ) 
{s±.}  j=l  i=l  10  13 


subject  to 


I 

I 

i=l 


*  T 


for  some  investment  target  T. 

22.  The  first  problem  will  be  considered  here.  When 

solving,  only  minor  changes  in  computing  give  the  solution 
of  this  second  problem. 


MARGINAL  ANALYSIS 


23.  Consider  a  fixed  item  i.  As  a  function  of  s^, 

log  A(s^q  ,  s^j)  is  a  concave  function,  where  s^Q  is  held 
fixed  (see  Annex  D) .  Thus  for  each  depot  stock  level  s^g 
of  a  fixed  item  i  one  can  obtain  a  maximal  value  of 

J  J 

Z  log  A(s.n  ,  s. .)  for  a  given  base  stock  S  =  E  s. .. 

j=l  1U  13  j=l  13 

In  fact  as  we  build  up  to  S  we  add  one  unit  at  a  time  giving 
it  to  that  base  which  yields  the  largest  increase  in  avail¬ 
ability.  Concavity  means  that  at  any  base  marginal  increases 
in  availability  for  each  new  addition  is  continually  decreas¬ 
ing.  Thus  if  a  unit  is  placed  at  a  base  which  gives  the 
largest  increase  in  availability,  all  subsequent  additions 
to  any  base  will  produce  lower  increases  in  availability. 

This  process  of  allocating  units  to  bases  one  at  a  time  does 
so  in  such  a  way  that  the  largest  marginal  increases  in 

E  log  A(si0  ,  Sj,j)  are  produced  first  and  when  S  is 
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attained  we  will  have  the  largest  possible  value  for 
J 

E  log  A  ( s  .  «  ,  s. .)  with  a  total  base  stock  S  and  depot 
j=l  1U  13 

stock  siQ  .  Call  this  value  SA(s^Q  ,  S) .  Now  for  each 
combination  of  base  stock  S  and  depot  stock  s^q  summing  up 
to  a  total  stock  m  (m  =  +  S)  pick  that  combination 

yielding  the  maximum  value  of  SA(s^q  ,  S)  and  call  this  value 
0^(m).  This  then  yields  the  best  possible  distribution  of 

m  units  of  item  i. 

24.  Now  we  must  perform  a  multi-item  distribution  in 

the  best  way  possible.  If  0 . (m)  were  concave  for  each  m 
then  the  same  arguments  as  above  would  apply  with  the 
exception  that  costs  must  be  accounted  for.  Thus  as  we  build 
up  the  total  stock  of  all  items,  at  each  step  the  next 
investment  is  allocated  to  the  item  that  produces  the  largest 
Increase  in  availability  per  dollar,  that  is  whose  value  of 


0i(m  +  1)  -  e^m) 


is  the  largest.  In  general  6^(m)  is  not  concave.  It  is, 
however,  a  montonically  increasing  function,  as  the  avail¬ 
ability  resulting  from  the  best  allocation  of  m  units  plus 

s  t 

an  arbi .  ary  allocation  of  the  m  +  1  unit  cannot  be  less 
than  tl .  xvailability  resulting  from  the  best  allocation  for 
m  unit.  It  is  also  bounded  from  above  as  A(s^q  , 
being  a  probability,  log  A(si0  ,  s^)  v<  0  so  that  e^m)  *  0. 
We  therefore  modify  0^(m)  by  defining  a  new  function  0|(m) 
so  as  to  take  values  on  the  boundary  of  the  convex  hull  of 
0^(m)  (smallest  convex  set  containing  the  graph  of  0^(m), 

see  Figure  2)  such  that  0£{m)  £  0^(m).  The  effect  of  this 
is  to  raise  only  those  values  at  which  0^(m)  is  not  concave 
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to  obtain  concavity.  So  now  we  would  allocate  the  next 
investment  to  the  item  i  with  the  largest  value  of 

0|(m  +  1)  -  8£<m) 


Allocation  terminates  when  the  total  availability  just 
exceeds  the  target  a. 

25.  Note  that  for  a  given  item  i,  we  don't  terminate 

allocation  at  a  level  m  such  that  0^(m)  <  0^(m) ,  as  the 
true  availability  of  item  i  would  be  less  than  the  value 
upon  which  a  decision  to  take  the  m*"*1  unit  was  made.  In 
fact  if  0i(m)  <  6|(m)  then  the  points  on  the  graph  of  0|  at 
the  values  m-1,  m  and  m+1  are  collinear  by  construction  of 
This  means  that 

0|(m+l)  -  0j(m)  =  0j(m)  -  0|(m-l) 

Hence 

0|  (m+1)  -  0j(m)  =  0!(m)  -  0£(m-l) 


Thus  the  next  best  investment  is  again  a  unit  of  item  i 
giving  the  same  increase  in  0£  per  unit  dollar.  Successive 
increments  in  item  i  are  then  continued  until  a  value  of  m 
is  reached  such  that  0|(m)  =  0^(m)  (which  must  occur  as  0^  is 
bounded  above).  For  example  in  Figure  2,  if  the  second 
unit  of  item  i  was  the  last  to  be  invested  in  then  a  third 
and  a  fourth  are  bought.  At  this  point  0|(4)  =  8|(4).  The 
point  here  is  that  investment  decisions  made  on  the  artificial 
values  of  0|  is  acceptable  as  they  point  to  future  large 
increases  of  0^  which  otherwise  would  be  missed.  Again  in 
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Figure  2  the  actual  increase  in  6,  going  from  1  to  2  is  so 
low  that  the  comparison  to  increases  incurred  by  other  items 
for  future  investments  may  never  permit  purchase  of  the 
second  unit  and  hence  we  may  never  benefit  from  the  large 
increase  obtained  in  going  from  3  to  4  units  of  item  i. 

26.  There  is  the  potential  problem  of  overshooting 
the  target  under  this  method  but  in  practice  this  does  not 
appear  to  be  a  serious  problem.  More  serious  is  the  fact 
that  with  a  large  number  of  items  (>  300)  the  computations 
of  the  multi-item  marginal  analysis  are  extremely  time 
consuming  on  a  computer.  The  next  section  describes  a 
computationally  more  efficient  method. 

THE  GENERALIZED  LA GRAN GIAN  METHOD 

27.  The  general  form  of  the  original  problem  is  the 

minimization  of  a  cost  function  H(x)  over  a  set  S  of 
strategies  given  constraints  c..,  j=l,  — ,  n  on  the  production 
Cj(x),  j=l,  - ,  n  of  n  products;  i.e. 

min  H (x) 
x  e  S 

subject  to 

Cj  (x)  *  Cj  j=l , - ,n 

The  following  theorem  by  Everett  (10)  will  simplify  the 
problem  by  changing  it  to  an  non-constrained  problem. 

Theorem  1 «  Let  t.,  j=l, — ,n  be  non-negative  real  numbers. 

n 

If  x*eS  minimizes  the  function  H(x)  -  £  r.  c.  (x)  over  all 

*  j»l  3  3 

xeS,  then  x  minimizes  H(x)  over  the  subset 
S*  *  (x|Cj(x)  >  Cj(x*)for  all}  cs 
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Proof ;  Let  x*  be  as  stated.  Then 

n  n 

H(x*)  -  Z  t.  c  .  (x*)  4  H (x)  -  Z  t.  c . (x)  for  all  xeS 
j=l  33  j-1  3  j 

n 

H  (x*)  *  H  (x)  +  Z  t  .  { c  •  (x* )  -  c.  (x)}  for  all  xeS 
j=l  3  3  3 

In  particular  this  last  inequality  is  true  on  S*CS.  However 
on  S* 

Z  t  .  {c  .  (x*)  -  c  .  (x)  }  $  0 

j  J  3  3 

as  >  0  for  all  y  But  this  means  that 

H(x*)  4  H (x)  on  S*  □ 

Notice  that  when  solving  the  unconstrained  problem  we  solve 
a  constrained  problem  which  is  more  restricted  than  the 
original  problem  and  of  course  it  is  not  known  in  advance. 

There  is  no  guarantee  that  x*  exists  or  is  unique  for  the 
unconstrained  problem.  If  one  such  x*  does  exist  it  gives 
the  minimum  cost  possible  without  producing  less  (availability 
in  our  case)  than  x*  does.  However,  running  through  a 
spectrum  of  values  for  ij's  gives  answers  under  a  series  of 
different  constraints  and  permits  the  values  of  (c.(x*)}  to 
approach  those  of  {c^}  through  trial  and  error.  Such 
evaluations  under  a  series  of  different  constraints  is 
common  procedure  in  operations  research. 

28.  Theorem  1  permits  us  to  reformulate  the  problem: 

for  a  given  t,  find  a  stock  strategy  optimizing 

I  J  II 

min  {  Z  (c.  Z  s.  .)}  -  x  Z  {  Z  log  A(s.n  ,  s4^)} 

{8ij}  i=l  1  j-0  ij  j-1  i-1  i0 


m 
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Here  we  find  the  second  important  feature  of  theorem  1  for 
the  form  of  our  problem,  namely  the  function  to  be  minimized 
can  be  separated  into  independent  functions ,  one  for  each 
item;  i.e. 

I  J  J 

min  E  {c.  E  s .  .  -  t  E  log  A(s..  ,  s.  .)} 

{s..}  i=l  1  j-0  13  j=l  10  13 


and  the  minimum  can  be  obtained  by  finding  the  minimum  for 
each  item  i  of 


F . 
l 


J 

-  t  E  log  A(s .  n  ,  s.  .) 
j=l  13 


Let  t  be  a  fixed  positive  number,  called  the  Lagrange  multiplier, 

J 

and  i  a  fixed  item.  For  each  m  «  E  s,  ., 

j=0  13 

Fi(m)  «  ci  m  -  t  6i(m) 

gives  the  minimum  value  of  Fi  for  that  value  m  of  total  stock. 

The  next  theorem  tells  us  that  there  is  a  unique  minimum  of 

F.  as  a  function  of  m  and  conditions  under  which  to  find  it. 
i 

Theorem  2  There  exists  a  unique  optimum  m  of 

min  F^(m) 
m 

satisfying 

ci  -  t  {6|(m)  -  e^m-l)}  <  0  4  C£  -  t  {6j!(m+l)“  6£(m)) 
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where  0!  is  the  concave  extention  of  0.  described  earlier, 
x  i 

For  this  m,  Q^(m)  =  0^(m). 

The  proof  of  this  result  is  found  in  Annex  C  where  -tl).  is 
taken  as  0  a  monotonically  decreasing  function  bounded 
below  by  0  (as  x  >  0  and  9^  is  monotone  increasing  and 
bounded  above  by  0) .  Thus  each  item  is  treated  individually 
(avoiding  the  multi-item  approach)  and  its  total  stock  m 
is  determined  by  the  conditions 

9 1  (in+1)  -  9[(m)  ±  9j(m)  -  9|(m-l) 


which  are  the  above  conditions  restated.  In  words,  we  stop 
allocating  stock  of  item  i  when  the  increase  in  availability 
per  unit  dollars  drops  below  — .  These  conditions  occur  at 
an  m  where  0j  (in)  actually  evaluates  the  total  availability 
0i  (in)  . 


29.  The  major  difficulty  with  the  Lagrange  method  is 

that  the  t  is  unknown.  One  could  display  several  cost- 
effectiveness  alternatives  with  different  values  x  or  if 
an  availability  target  is  known,  judicious  trial  and  error 
runs  can  be  carried  out  as  described  earlier. 


IV  -  COMPUTATIONS 


30.  In  this  section  we  describe  the  computational 

steps  required  for  the  solution  of  the  optimization  problem 
described  in  section  III.  If  the  vector  (s’Jj’g  ,  s*^,  ...,  s *Jj\) 

is  the  best  allocation  of  m  units  of  item  i  as  described  in 
section  III  then  the  availability  of  item  i  was  given  by 


0. 

1 


.in 


J 

(m)  =  Z  log  A(s“  ,  s‘"  ) 
j=l  13 


The  availability  A(s™g  ,  s™.)  is  determined  by 


A  (s) 


s 

=  Z 
x=o 


P(x| AT) 


at  each  base  j  and  for  each  item  i  as  described  in  section 
II  para.  17  where  s  is  replaced  by  s”j  and  T  by  Tj  (siQ),  the 
mean  resupply  time  at  base  j  when  the  depot  stock  is  siQ. 

The  formal  steps  are  described  next.  Of  course  any  expression 
in  the  following  involving  P(x|Xt)  is  modified  using  the 
Bayesian  procedure  described  in  paragraph  19  section  II, 
when  observed  demand  data  is  available. 

Step  1.  For  each  item  i  and  depot  stock  level  calculate 

the  resupply  time  t^(s^q).  Since  a  fraction  r^  of  the  demands 
are  processed  at  base  j ,  if  the  average  base  repair  time  is 
Aj  then  the  contribution  to  resupply  time  at  base  j  would  be 
r^-Aj.  If  an  item  leaves  base  j  to  be  repaired  at  the  depot 
the  time  required  would  be  the  order  and  shipping  time  0^  if 
there  was  a  infinite  amount  of  depot  spare  stock,  but  would 
be  0^  +  D,  where  D  is  the  average  depot  repair  time  (D 
actually  consists  of  the  average  time  for  an  item  to  be 
shipped  from  base  to  depot  and  to  be  repaired  at  depot) ,  if 
that  spare  stock  is  zero.  Thus  the  actual  resupply  time  is 


where  the  function  6{s^q)  represents  the  ratio  of  expected 
number  of  units  incurring  delays  during  period  D  to  the 
expected  total  demand  during  period  D.  As  demands  come  in 
from  all  bases  the  customer  arrival  rate  is 

J 

A  =  I  (1-r.)  A .  a  . , 
j=l  J  J  J 

where  a^  is  the  change  in  program  factor  from  data  period 

to  prediction  period  at  base  j.  This  follows  from  Theorem 

3  of  Annex  A  and  the  fact  that  the  customer  arrival  rate 

from  base  j  is  (1-r.)  A.  a..  Thus  with  a  stock  level  of  s.n 

3  3  3 

for  item  i  at  the  depot,  the  expected  number  of  delays  in 
time  D  is  the  expected  backorder  (see  para.  16) 


B(si0|AD)  =  Z  (x“si0)  P(*1XD) 

x~si0 


and  the  expected  total  demand  in  time  D  is  XDf. 


6Ui0> 


B(si0|AD) 

ADf 


Thus 


As  the  fraction  of  items  going  to  the  depot  for  repair  is  1-r 
the  contribution  to  the  resupply  time  for  these  items  is 


(1-r.)  (0.  +  6(si0)  D) 

Finally  we  compute 


T.(s.J  *=  r .  A,  +  (1-r.)  <0.  +  5  (s . «)  D) 
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Step  2.  For  each  depot  stock  level  s ^  compute  the 
availability  for  all  stock  levels  s ^  at  base  j.  Thus  we 
calculate  for  various  levels  of  s_ 


A(si0 


V  = 


A  ( s  .  .|X.  T  . 
ID1  D  J 


(si0)  > 


In  this  way  we  obtain  the  following  table  of  availability 

for  each  i  and  value  s . n 

lO 


This  table  is  calculated  up  to  the  maximum  allowed  level  of 
stock  MBST  for  item  i  at  base  j. 


Step  3. 


iO 


S)  ■  max  I  log  A(s 
Esij"S  ^ 


iO 


sij) 


Compute  SA(s 


9 
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With  s^q  fixed,  then  according  to  section  III  paragraph  23, 

units  of  an  item  i  can  be  distributed  one  at  a  time  where 
the  next  unit  is  added  to  that  base  where  the  largest 
increase  in  availability  will  occur.  Thus  we  can  form  the 
following  table  of  best  availability  possible  for  a  stock 
level  S  of  item  i  when  the  depot  has  siQ  stock  of  item  i. 
The  calculations  are  done  up  to  the  maximum  system  stock 
level  MSTAL  allocatable  to  bases  for  item  i  and  for  each 
depot  stock  level  up  to  its  maximum  MDST. 


Step  4.  Compute  0^(m) 


max  SA(s.ft  ,  S) 
si0+S«m 


For  these  computations  we  look  along  the  diagonals  of  the 
table  in  Step  3  defined  by 


si0  +  S  *  constant  *  m 

and  choose  the  maximum  value  in  the  diagonal.  This  is  the 
best  allocation  of  m  units  of  item  i. 
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Step  5»  Marginal  Analysis  for  the  multi-item  problem. 

As  0^(m)  may  not  be  concave  the  concave  extension  0£  of  0^ 

must  be  performed  first.  Next  A  0^(m)  =  0^(m)  -  0.?(m-l)  is 

calculated.  At  each  step,  the  next  investment  is  allocated 

to  the  item  that  produces  the  largest  value  of  -  where 

c.  is  the  unit  cost  of  item  i.  Allocation  terminates  when 
1 

an  investment  target  or  a  system  availability  target  is  just 
exceeded.  If  an  optimal  allocation  of  stock  using  the 
Lagrange  multiplier  method  is  required  then  an  item  by  item 
allocation  is  performed,  where  allocation  terminates  for  an 
item  i  when  a  total  stock  level  m  is  reached  such  that 

c.  c . 

. -  -  J: _  <  T  «  1 

A  eT{m)  F ' 6]  (m+1) 

First  estimates  for  t  may  be  obtained,  for  instance,  by 
taking  a  small  item  sample,  solving  using  the  multi-item 
approach  and  if  m  is  the  optimal  stockage  policy  for  a 
sample  item  i  use 


A8|  (m) 


or  average  this  value  over  the  sample,  etc.  (Note  that 
budget  or  availability  target  values  would  have  to  be  down 
graded  proportionately  to  the  sample  size.) 
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31.  There  are  several  different  ways  to  use  this 

sequence  of  calculations.  If  we  want  an  evaluation  of  the 
performance  of  existing  stock  and  investment  cost, 
computations  can  stop  at  step  2.  If  a  redistribution  of 
existing  stock  to  improve  availability  is  desired  then 
steps  1  to  4  will  optimally  reallocate  total  assets  among 
the  bases  and  the  depot.  In  the  five  step  calculations 
minimum  stock  levels  were  taken  as  zero.  If  minimum  stock 
levels  of  an  item  at  base  or  system  level  are  specified 
then  computations  for  values  lower  than  these  levels  will 
be  omitted.  After  the  above  allocation  is  completed  the 
total  system  investment  and  system  availability  are  computed. 
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ANNEX  A 


FRSC1DINQ  PiQB  BUNK-NOT  HU<jED 


THE  COMPOUND  POISSON  PROCESS 


Here  we  derive  the  mean  and  variance  of  the  compound 
Poisson  distribution  and  show  that  the  logarithmic  Poisson 
process  has  the  probability  states  of  the  negative  binomial 
distribution.  To  this  end  the  concept  of  the  characteristic 
function  is  useful  in  establishing  relationships.  If  the 
random  variable  X  has  a  probability  distribution  function  g(x) 
then  the  characteristic  ijjg(u)  is  defined  by 

where  u  is  a  real  number,  i  =  /-I  and  E  is  the  expected  value. 

We  note  that  distinct  probability  distributions  have  distinct 
characteristic  functions  (see  [l2;XV.3,  Theorem  l]  ). 

Lemma  1 .  Let  X  and  Y  be  two  independently  distributed 
random  variables  with  probability  distribution  functions  f 
and  g  respectively.  Then 

<J>f*g(u)  =  1»f(u)  •  ^g(u) 


where  f*g  is  the  convolution  of  f  and  g  (i.e.  the  p.d.f.  of 
X  +  Y). 

Proof:  ^f*g(u)  =  E(elu(X+Y))  =  E(eluX  •  eiuY) 

=  E  (eluX)  •  E  (eluY)  =  i|»f(u)  •  (u) 


o 
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Theorem  1 .  Let  P(x|Xt)  be  the  p.d.f.  of  a  compound 

Poisson  distribution  with  demand  p.d.f.  f(w)(see  para  12).  Then 

♦  ptu)  =  eU  '♦f'"'-1' 

oo 

Proof:  ij/_(u)  =  E  (e1U  X(t))  =  E  P(x|Xt)  elux 

x=0 


=  E 
x=0 


I  mi  e-Xt  fy*(x) 

y=0  y ' 


1UX 


e"Xt  E  (Xt)y  E  fy*(x)  elux 
y=0  y!  x=0 


e-xt  z  iXt)_  ^  yMu) 

y=0  y! 


=  e  E 


y=0  yl 


by  lemma  1 


_  -Xt  Xt  <u) 
—  e  e  x 


_  gXt  (<Jjf(u)-l) 


o 
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Remark:  The  change  in  the  order  of  summation  in  the  proof 
is  valid  since  for  any  fixed  x  the  sum  over  the  y's  is  a 
finite  sum  (see  [2§ ;  theorem  8.3]). 

Theorem  2  For  P(x|At)  we  have 

m  =  At  E (W) 

o2  =  At  E(W2) 

Proof  From  the  definition  of  (u) 

00 

d  VU)  L_n  =  1  P(x  |  At)  ix  =  i  E  (X)  =  im 

du  u  u  x=0 

d2  ^p(u)  00  j  _ 

- , -  ln  =  E  P  (x  |  At)  (ix)  z  =  -  E  (X  ; 

du  u-u  x=0 


On  the  other  hand  from  the  expression  of  tf>p(u)  in  the 
previous  theorem 


d  <J>p(u) 
du 


u=0 


eAt  (*f (o)-l)  Xt  d^f(u) 

du 


u=0 


=  At  i  E (W) 


as  i|>f(o)  =  1 


r 


d  t|»p(u) 

~7~l - 

du 


(\t  d^f  (u)  (  _  XtO)j^(o)-l) 

u=0=V— —  u=0  • 


At  (4>^ (o) -1)  d  *f(u)  , 

!  f  At  ' u=0 

du 


=  (At  i  E(W) ) 2  +  At  (-E(W2) ) 
=  -At  [At  E(W)2+  E(W2) | 

Therefore 

im  =  At  i  E (W) 
m  =  At  E(W) 

and  E(X2)  *  At  [  At  E(W)2+  E(W2) ] 
so  that  a2  =  E(X2)  -  m2 
=  At  E(W2) 

a 


At  the  depot  compound  Poisson  demands  are  coming  in  from  all 
bases.  That  the  net  process  at  the  depot  is  again  compound 
Poisson  follows  from  the  following  theorem. 


Jlipl.il  I II I  !■ 


-  A5  - 


Theorem  3  .  A  sum  of  N  independent  compound  Poisson 

process  with  arrival  rates  A^  and  compounding  functions 

f . ,  i=l,  N  is  a  compound  Poisson  process  with  arrival 

1  N 

rate  A  =  £  A .  and  compounding  function  f  such  that 
i=l 

N  (u) 

i|;_(u)  =  E  — - - - 

r  i=l  * 

Proof:  Since  this  is  true  for  N=l,  by  induction  assume  it  is 

k+1 

true  for  N£k.  Consider  the  sum  X  =  Z  X.  of  k+1  independent 

i=l  1 


compound  Poisson  processes.  By  assumption  then  X’  =  £  X.  is  a 

i=l 

k 

compound  Poisson  process  with  arrival  rate  \  -  Z  A.  and 

i«l  1 

compounding  function  f  such  that 

*  .  <u) 

«Pf(u)  =  Z  _ 

i=l  A 


If  P  is  the  probability  distribution  for  X  =  X'  +  X^+^  then 
by  lemma  1  and  theorem  1 

(u)-l) 

*p(u)  -  •  e 


k+1 


=  e  *  *  +  Xk+i  *fk+1  "  W 


t(  E  Ai  ip  +  Ak  ip  -  A  -  A  ) 
i=l  1  ri  K 1  rk+l  K+1 


=  e 


k+1 

(  Z  A.)  t 
i=l 


/k+1  /k+1  \ 

[  E  A.  /  E  A. I  - 

\i«i  1  fi/i-i  V 


e 
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As  X  is  the  sum  of  two  compound  Poisson  processes  X'  and  X^+^, 
it  is  by  assumption  compound  Poisson  and  by  uniqueness  of 
representation  by  characteristic  functions,  its  arrival  rate 

k+1 

must  be  £  X.  and  the  characteristic  of  the  compounding 

1=1  / 

k+1  /k+1 

distribution  must  be  I  X,  ip-  /  Z  X.  .  This  shows  the 

i=i  1  *y  i=l  1 

statement  is  true  for  N  =  k  +  1.  ° 


Finally  it  is  shown  that  the  probability  distribution  P(x|Xt) 
of  a  logarithmic  Poisson  process  is  a  negative  binomial 
distribution.  The  negative  binomial  distribution  gives  the 
probability  of  the  rfc  success  at  the  (r+k)  trial,  namely 
if  s  is  the  probability  of  success  and  t  «  1-s  the  probability 
of  failure  then 


where  by  convention  ( x\  =  — ijl: (x-r+1)  is  extended  to 

W  r! 

permit  x  to  be  any  real  number.  The  name  of  the  distribution 
comes  from  the  fact  that 


Note  that  for  r=l  this  is  the  geometric  distribution. 
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By  uniqueness  of  representation  through  characteristic 
functions  it  suffices  to  show  that  the  characteristic 
functions  of  these  distributions  are  identical.  Recalling 
that  the  compounding  distribution  of  the  logarithmic  process 
is  given  by 


then 


f  (w)  = 


0 

'  1 
w  An  q 

JL _ 

in  q 


Z 

w=l 


w  =  0 

w  3  1,2,3,... 
q  =  p+1  >  1 


_i _  [-An<l-  |  e1U>] 

m  q 


X  [in  q  -  An  (q  -  PelU)] 

An  q 

An  (q  -  peiu) 

An  q 


so  that 


Vu> 


eXt  (u) -1 
ek(-An(q-peiu)  ) 


(Xt  «  k  tn  q) 


k 


ANNEX  B 


PROOF  OF  PALM'S  THEOREM 


Let  F ( t) 


^(f;)d£  be  the  distribution  function  for  4> 


Proposition  1.  The  probability  that  any  particular  order  placed 


in  the  interval  0  to  t  has  arrived  by  time  t  i 


ts 


FU)d£. 


Proof ;  First  we  consider  the  probability  that  a  demand  occurs 
between  5  and  5  +  d£  knowing  that  it  occurred  in  [o,t3-  This 
is  given  by 

Pr(l  demand  in  (k»€+dc])«  Pr(0  demands  in  [o,€]u[c+d€,t]  ) 

Pr  (1  demand  in  [0,t3) 

-Ad£  -A ( t-d£) 

=  Adge _ .  e _  =  d£ 

-At  t 

Ate 

Next  the  probability  that  a  unit  ordered  at  time  £  will  arrive 
by  time  t  is  F(t-£).  Thus  knowing  that  at  least  one  demand 
occurred  in  [0,t]  ,  the  probability  that  it  occurred  between 

£  and  d£  and  that  the  unit  ordered  arrived  by  time  t  is 

F(t-C)  d£ 

€ - 

Consequently  the  probability  that  any  particular  order  placed 
in  the  interval  0  to  t  will  arrive  by  time  t  is 

d5 

An  obvious  change  of  variable  yields  the  result  of  the 
proposition. 


Proposition  2.  The  expected  resupply  time  can  be  expressed 
as 


T  =  (l-F(U)d^ 


Proof ; 
T 


By  definition 


CdF(C) 


B2 


We  use  integration  by  parts  to  obtain 


/ 


b 

CdF(C)  =  bP(b) 

0 


F(£)d£ 


=  b(F(b)  - 


But  0  <  b (1  -  F (b ) ) 


which  tends  to  0  as  b— - «►“,  giving  the  desired  result. 


Proposition  3  The  steady  state  probabilities  of  x  units  in 
resupply  are  given  by  the  Poisson  process  with  rate  AT. 

Proof ;  Suppose  the  system  has  been  operating  for  time  t.  If 
during  this  time  n  demands  were  placed  the  number  x  of  units 
in  resupply  at  time  t  is 
x  =  n  -  arrivals 

Thus  n  -  x  would  have  arrived.  Under  these  circumstances  the 
probability  of  x  units  in  resupply,  can  be  given  by  the 
binomial  probability  of  n  -  x  successes  (arrivals)  in  n  trials 


with  p  = 


F  ( {; )  d£  as  the  probability  of  success  (proposition  1) 


and  x  failures  with  probability  of  failure  q 


FU))d£  . 


Thus  the  probability  of  x  units  in  resupply  (failure)  from  n 
demancs  (trials)  is  B(x:n,q)  . 

As  the  probability  of  getting  n  demands  is  P(n|At)  in  time  t,  the 
probability  that  x  units  in  resupply  at  time  t  is 


Q  (x  1 1 )  =  Z  P(njAt)  B(x;n,q) 
n=x 


i. 
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«  l  (At)ne~Xt  .  n! 

'  xT  (n-x) : 


x  n-x 

q  p 


n=x  n; 

x  _x  -At 


T.  (At)n~X  pn_X 


=  (At)  q  e  _ 

xT  n=x  (n-x)  I 

to 

=  (Atci)X  e'At  £  (Lt£)n"X 
xl  n=x  (n-x) I 

, . ^  .x  -At  Atp 
=  (Atq)  e  e  ^ 


-  (Xta)x  e~xtq 

X.' 

=  P  (x I Atq) 


By  proposition  2 


,  tq  =1^  (1- 

J  o 


F(£))d£  T  as  t  — ►  00  so  that 


lim  Q (x | t )  =  P (x | AT)  showing  that  the  limiting  state 

t  — *00 

probabilities  of  x  units  in  resupply  are  Poisson,  depending 
only  on  the  mean  resupply  time  T. 


FPfr-TOl  uf»  PaOS  BLANK-NO!  FILLED 

ANNEX  C 


PROOF  OF  THEOREM  2  -  SECTION  III 


Here  we  prove  the  following: 

Theorem .  Let  0  be  a  real  valued  monctonical ly  decreasing  function 
defined  on  the  set  of  non  negative  integers  and  bounded  from 
below.  Let  0'  be  defined  on  the  non  negative  integers  taking 
values  on  the  boundary  of  the  convex  hull  of  6  such  that 
0 1 (m)  <  0 (m) .  if  c  >  0  there  exists  a  unique  minimum  m  of 

(1 )  min  cm  +  0 (m) 

m 

satisfying 

(2)  c+0 '  (m)  -  0'(m-l)  <  0  N<  c+0'(m+l)  -  0 '  (m) 

In  this  case 

0  '  (in)  =  0  (in) 

Proof:  0'  is  decreasing  as  any  two  consecutive  points  of  the 

graph  of  0'  are  on  a  line  joining  two  points  of  the  graph  of 
0  which  must  have  a  non-positive  slope.  Also  0'  is  bounded 
from  below.  In  fact  as  6  is  bounded  from  below,  by  say  b, 
then  the  graph  of  0,  G(9),  is  contained  in  B={ (x,y)  |y>b}  which 
is  a  convex  set  (see  Figure  1).  Thus  G(0')  C  {convex  hull  of 

G(0)  }CB  and  so  0'(m)  is  bounded  below  by  b  as  well.  Finally 

6'  is  clearly  convex  by  definition.  Thus  0 ' (m)  converges  and 
by  convexity  of  0',  4 (m)  =  0 ' (m)  -  0'(m+l)  is  monotonically 
decreasing.  If  we  define  O'  ( —  1 )  =  “>  then  A(-l)  =  «*.  So  by 

monotonicity  of  A(m)  there  is  a  unique  m,  call  it  m  >  0,  such 

that 

c  <  A(m-l) 

and 


c  >  A (m) 


m  =  1,  2,  3  are  all  minima  for  cm  +  6 ' (m) 

but  m  is  chosen  as  the  one  with  the  smallest  value 


C3 


or 

c  -  A(m-l)  <  0  4  c  -  A  (m) 

which  is  the  condition  (2)  stated  in  the  theorem.  Also  as 
cm+0 '  (m)  is  convex  m  is  a  a  minimum  of  cm+0 '  (m)  as  condition  (2) 
can  be  written  as 

c (m-1)  +  0 ' (m-1 )  >  cm  +  0 1 (in) 

cm  +  6 '  (m)  <  c(m+l)  +  0'(in+l) 

(Heuristically ,  all  minima  of  c  m+  1 (m)  must  be  adjacent  and 
condition  (2)  picks  out  the  one  with  the  smallest  value  -  see 
Figure  2.)  Next  we  show  that  6  (in)  =  0  '  (in)  .  If  0  (in)  ?  0  ’  (in) 
then  by  definition,  the  value  of  8'  at  in  has  been  lowered 
making  the  points  of  G(0')  at  m-1,  in  and  in+1  collinear.  That 
is 

0 1  (in)  -  0 1  (in-1)  =  0'(m+l)  -  0'(in) 

But  this  contradicts  condition  (2) .  Finally  as 

cm  +  0 (m)  =  cm  +  0 ' (m)  <  cm  +  0 ’ (m)  <  cm  +  0 (m)  Vm 

cm+0 (m)  has  a  minimum  at  m 


ANNEX  D 


PROOF  OF  THE  CONCAVITY  OF  LOG  A(s) 


The  availability  A(s)  used  in  METRIC  2  is  defined  by 
s 

A  (s )  =  7  P  (x  |  A  t ) 

x=0 

(see  paragraph  16/  section  II),  where  P(x]At)  is  the  logarithmic 
Poisson  distribution  (para  13,  section  II)  given  recursively  by 

P  (0  |  A t)  =  l/qk,  P  (x+1 1  At)  =  P  (X|  xt)  x=  1 , 2  ,  .  .  . 

and  where  s  is  the  stock  level  at  a  base.  For  brevity  we  will 

denote  P(x|At)  simply  by  P  (x)  and  by  C(x).  In  METRIC  2, 

■Log  A (s )  is  used  in  the  marginal  analysis  process  instead  of 
the  Rand's  expected  backorder  function.  Convexity  must  be 
assured  for  use  in  the  marginal  analysis  process.  The  expected 
backorder  function  is  convex  from  its  definition.  This  must  be 
verified  for  -  Log  A(s).  We  will  show  here  that  Log  A(s)  is 
concave . 


Before  going  on  to  the  proofs,  we  note  the  following, 
P(x+1)  =  C  (x)  P(x) 

and  C  (x)  -  C(x-l)  =  (x^^x  •  Thus  if  k>l,  C(x)  is  a 

monotonically  decreasing  sequency  converging  to  =  i-i. 
where  the  variance  to  mean  ratio  q  is  greater  than  1.  If  k>l, 
C(x)  is  a  monotonically  increasing  sequence  converging  to  1-^- 
and  if  k-l  C(x)  =  1-i  for  all 


D2 


The  next  two  theorems  give  the  proof  of  the  concavity 
of  Log  A (s) 

Theorem  1  (C(s)-l)  A(s)  <  P(s+1)  for  all  s  >  0. 

Proof :  We  look  at  two  cases. 

1)  If  k<l  then  as  C(s)  increases  to  1~-  we  have 

C(s)<l-—  so  that  C(s)-l<~i-<0  for  all  s.  Therefore, 
q  <3 

tC(s)-l)  A (s)  <  0  <  P(s+1)  for  all  s. 


2)  If  k>l  then  we  do  a  proof  by  induction  on  s. 

For  s*0 

(C(O)-DA(fl)  =  (C(0)-1)  P<0) 

«  C{0)  P(0)  -  P(0) 

=  P(l)  -  P  (0)  <  P(l) 

and  the  hypothesis  is  verified  at  s=o.  Assume 
it  is  true  for  s-1,  i.e. 

(C(s-l)-l)  A  (s-1 )  <  P{s) 

Since  C(s)  is  a  decreasing  sequence  then 
C(s)  <  C(s-l).  The  above  inequality  becomes 

(C(s)-l)  A(s-l)  <  P(s) 

P(s+1)  -  P(s)  +  (C(s)-l)  A(s-l)  <  P(s+1) 

(C(s)-l)  P(s)  +  (C(s)-l)  A (s-1)  <  P(s+1) 

(C(s)-l)  (P(s)  +  A(s-l) )  <  P(s-tl) 

(C(s)-l)  A(s)  <  P(s+1) 

which  is  the  hypothesis  for  s  a 


Theorem  2  Log  A(s)  is  concave. 

Proof:  To  show  that  Log  A(s)  is  concave  we  must  show 

A2  Log  A (s )  <  0,  i.e.  Log  A(s+2)-2  Log  A(s+l)+Log  A(s)<0 
Ffom  theorem  1  at  s+1  we  have 
(C (s+l)-l)  A(s+1)  <  P  (s+2) 
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Multiplying  both  sides  by  P(s+1)  we  obtain 
(P(s+1)  C(s+1)  -  P  (s+1) )  A  (s+1 )  <  P(s+1)  P  (s+2) 

(P  (s+2 )  -  P(s+1))  A  (s+1 )  -  P(s+1)  P  (s+2 )  <  0 

(A (s+1) )  2+  (P (s+2 ) -P  (s+1) )  A ( s+1 ) -P (s+l)P (s+2) < (A (s+1) ) 2 

(A(s+1)+P (s+2) )  (A (s+1) -P ( s+1) ) < (A (s+1 )  2 

A  (s+2 )  A  (s )  <  (A  ( S+1 )  )  2 

Log  A (s+2)  +  Log  A(s)  <  2  Log  A (s+1) 

Log  A (s+2 )  -  2  Log  A(s+1)  +  Log  A(s)  <  0 
which  is  the  required  result.  D 
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